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1 Introduction 



The main goal of this paper is to construct the integrable model which corresponds to 
the Af = 2 SUSY SU(N) Yang-Mills theory with a hypermultiplet in the antisymmetric 
representation. The 1994 work of Seiberg and Witten [|lj had shown that the Wilson effective 
action of M = 2 SUSY Yang-Mills theory is determined by a fibration of spectral curves T 
equipped with a meromorphic one-form dX, now known as the Seiberg- Witten differential. 
It was soon recognized afterwards [@, [3], |J that this set-up is indicative of an underlying 
integrable model, with the vacuum moduli of the Yang-Mills theory corresponding to the 
action variables of the integrable model. In fact, in the special case of hyperelliptic curves, a 
similar set-up for the construction of action variables as periods of a meromorphic differential 
had been introduced in ||. This unexpected relation between TV = 2 Yang-Mills theories on 
one hand and integrable models has proven to be very beneficial for both sides. The Seiberg- 
Witten differential has led to a universal symplectic form for soliton equations in the Lax or 
Zakharov-Shabat representation [g, [7j. The connection with integrable models has helped 
solve the SU(N) Yang-Mills theory with a hypermultiplet in the adjoint representation [f|. |J. 
as well as pure Yang-Mills theories with arbitrary simple gauge groups G || . Conversely, the 
connection with Yang-Mills theories has led to new integrable models, such as the twisted 
Calogero-Moser systems associated with Yang-Mills theories with non-simply laced gauge 
group and matter in the adjoint representation ||, and the elliptic analog of the Toda lattice 

E3- D 

Despite all these successes, we still do not know at this moment how to identify or 
construct the correct integrable model corresponding to a given Yang-Mills theory. This is a 
serious drawback, since the integrable model can be instrumental in investigating key physical 
issues such as duality, the renormalization group, or instanton corrections |13|, [TJ], [D|. At the 
same time, the list of spectral curves continues to grow, thanks in particular to methods from 
M theory |L6], [T7| and geometric engineering ]TB) . It seems urgent to develop methods which 
can identify the correct integrable model from a given spectral curve and Seiberg- Witten 
differential. 

In the case of interest in this paper, namely the SU(N) gauge theory with antisymmetric 
matter, the Seiberg- Witten differential and spectral curve had been found by Landsteiner 
and Lopez |DJ using branes and M theory. The Seiberg- Witten differential dX is given by 



dX = x—. (1.1) 

y 

The spectral curve is of the form 

y 3 _ (3A Af+2 + x 2 Y^ u . x i) y 2 + ( 3A ^+2 + x 2 <£ { _y ^ ^ +2 y _ A 3(Af + 2) = Q> (L2) 
i=0 i=0 

where A is a renormalization scale. For the SU(N) gauge theories, one restricts to un = 1, 
un-i = 0, so that the moduli dimension is N—l, which is the rank of the gauge group SU(N). 



1 We refer to O, [T^] for more complete lists of references. 



2 



The Landsteiner-Lopez curve ( |1.2| ) and differential ( |1.1| ) have been studied extensively by 



Ennes, Naculich, Rhedin, and Schnitzer [20]. In particular, they have verified that the curve 
and differential do reproduce the correct perturbative behavior of the prepotential predicted 
by asymptotic freedom. The problem which we wish to address here is the one of finding 
a dynamical system which is integrable in the sense that it admits a Lax pair, and which 
corresponds to the Landsteiner-Lopez curve and Seiberg-Witten differential ( |1 . 1| ) in the sense 



that its spectral curve is of the form (|1.2|) , and its action variables are the periods of dX along 
N — 1 suitable cycles on T. 

We have succeeded in constructing two integrable spin chain models, whose spectral 
curves are given exactly by the Landsteiner-Lopez curves. However, the action variables of 
the desired integrable model must be given by dX = x-^-, and here the two models differ 
significantly. For one model, referred to as the odd divisor spin model, the 2-form resulting 
from dX vanishes identically. For the other, referred to as the even divisor spin model, 
the Hamiltonian reduction of the 2-form resulting from dX to the moduli space of vacua 
{un = 1, un-i = 0} is non-degenerate, and the reduced system is indeed Hamiltonian 
with respect to this symplectic form, with Hamiltonian H = u^_ 2 - Thus the latter model is 
the integrable system we are looking for. 

Our main result is as follows^]. Let q n ,p n be 3-dimensional vectors which are N + 2 
periodic, i.e. p n+ N+2 = Pn, q n +N+2 = <?n, and satisfy the constraints 



Pn = goP-n-1 

where go is the diagonal matrix 
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Consider the dynamical system 

Pn+l Pn-1 



Pr 



+ 



T 1 T 

Pn+lQn Pn-lQn 



f^nPni Qr 



Q.n+1 



gn-1 

PlVn+i Plq n -i 



t^n Pn ■ 



for some scalar functions /i„(t). The system is invariant under the gauge group G 
by the following gauge transformations 

Pn * X n p n , Q n > X n Q n , 

Pn -> W T p n , q n -> W~ l q n , 

Here W is a 3 x 3 matrix which commutes with g , Wgo = g$W . Define the 3x3 
L(x) and M(x) by 



(1.6) 
generated 

(1.7) 
(1.8) 
matrices 



N+l 

L ( X ) = II ( 1 + X 1riP 
n=0 



'n ) j 



M(x) 



X 



( Qn+ipI _ QoPn+i 
\PoQn+i Pn+iQo, 



1.9) 



2 The notation is explained in greater detail in §3 and §5. 
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Main Theorem. • The dynamical system ( |1.6| ) is equivalent to the Lax equation 

L(x) — [M(x), L{x)]\ (1.10) 

• The spectral curves Y = {R(x, y) = det(yl — L(x)) = 0} are invariant under the flow (|1.6|) , 
and are exactly the curves of the Landsteiner-Lopez form (with A N+2 normalized to 1); 

• There is a natural map (q n ,Pn) — > (I\ D) from the space of all spin chains satisfying the 
constraints ([0|[L4]) to the space of pairs (T,D), where T is a Landsteiner-Lopez curve, and 
D = {zi, ■ ■ ■ , z 2 n+i\ is a divisor whose class [D] = [D a ] is symmetric under the involution 

a: (x,y) = z^z° = (-x,y- 1 ). (1.11) 

For a given (q n ,Pn), D is the set of poles of the Bloch function ip , L(x)ip = yip (x); 

• Let A4q be the space of pairs {T, [D]}, where T is a Landsteiner-Lopez curve with un = 1, 
un-i = 0, and [D] is a divisor class which is symmetric under the involution a. Then the 
space Ado has dimension 2(N — 1). The map (q n ,p n ) — > (r, D) descends to a map between 
the two spaces 

{{qn,Pn)}/G^M , (1.12) 

where on the left hand side, we have factored out the gauge group G from the space of periodic 
spin chains satisfying the constraints (|1.3J1-4|) . At a generic curve T and a divisor [D] in 
general position, the map ( |1.12| ) is a local isomorphism. 



Let the action variables ai and the angle variables (pi be defined on the space A4 by 

2N+1 z . 

dX, 0i = E / <k>i (1-13) 
A i i=i J 

where {A}i<i<7v-i and {dui}i<i<N_i, are respectively a basis for the even cycles and a basis 
for the even holomorphic differentials on T. Then 

N-l 

uj = ^2§ ai A 5<pi (1.14) 

i=l 

defines a symplectic form on the 2(N — 1) -dimensional space A4q, 

• The dynamical system ( |1.6| ) is Hamiltonian with respect the symplectic form ( |1.14| ) . The 
Hamiltonian is H = Un- 2 - 

In terms of the {q n ,Pn) dynamical variables, the Hamiltonian can be expressed under the 
form 

H _ U N _ 2 U 2 N _ t _ ™ (Pngn-3) (Pngn-2) 2 ^ 

u N 2u 2 N ^ (plq n -i){Pn-iQn-2){pl-2Qn-3) 2(p^g n _i) 2 (pl-tfn-2) 2 ' 

u 2 

where we have used the constraint u N = 1, = to write H as H = — ^j- 1 



We would like to note the similarity of the Lax matrix L in (1.9) to the 2x2 Lax matrix 



used in |21[ for the integration of a quasi-classical approximation to a system of reggeons in 
QCD. 
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A key tool in our analysis is the construction of |], |7|, which shows that symplectic 
forms constructed in terms of Seiberg-Witten differentials can also be constructed directly in 
terms of the Lax representation of integrable models. The latter are given by the following 
universal formula || [7| 

u = - Res Pa < i)* n+l 5L n (x) A 5ip n > dx (1.16) 

where ip n and ipn+i are the Bloch and dual Bloch functions of the system, and P a are marked 
punctures on the spectral curve Y. In the present the 3 points on Y above x = oo. 

Finally, we note that the odd divisor spin model (which we describe in §3.1 and §6) may be 
of independent interest. Although the symplectic form associated to the Seiberg-differential 
x^- is degenerate in this case, the model does admits a Hamiltonian structure with non- 

degenerate symplectic form, but one which is associated rather with the form dA(i) = lay— . 
As suggested in ||19|| , the form Iny^- is also indicative of supersymmetric Yang-Mills theories, 



but in 5 or 6 dimensions with M = 1 supersymmetry. 



2 Geometry of the Landsteiner-Lopez curve 

We begin by identifying the geometric features of the generic Landsteiner-Lopez curve which 
will play an important role in the sequel. Fixing the normalization A N+2 = 1, we can write 

T : R(x, y)=y 3 - f(x)y 2 + f(-x)y -1 = (2.1) 

where f(x) is a polynomial of the form 

N 

f(x)=3 + x 2 P N (x), P N (x) = J2^ 1 (2-2) 

The parameters Uo, • • • , Wjv are the moduli of the Landsteiner-Lopez curve. 

• The Landsteiner-Lopez curve T is a three-fold covering of the complex plane in the x 
variable. It is invariant under the involution a defined in ( |1 . 1 1|) . The important points on Y 
are the singular points, the points above x = oo, and the branch points. We discuss now all 
these points in turn. 

• The singular points are the points where 

d x R(x,y) = d y R(x,y) = (2.3) 

The generic Landsteiner-Lopez curve has exactly one singular point, namely (x, y) = (0, 1). 
At this point, the equation ( |2.1| ) has a triple root, and all three sheets of the curve intersect. 
For generic values of the moduli Ui, all three solutions y of R(x,y) = can be expressed as 
power series in x in a neighborhood of x = 

oo 

y{x) = l + Y j y^ 1 (2-4) 

i=l 
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In fact, we can substitute ( |2.4| ) into ( |2.1| ) to find recursively all coefficients y iy with the first 
coefficient y x a solution of 

y\ ~ u yi + 2ui = 0. (2.5) 

For generic Uq, U\, this equation does admit three distinct solutions for yi, which lead in turn 
to the three distinct solutions. These three distinct solutions provide effectively a smooth 
resolution of the curve T, where the crossing point y = 1 above x = has been separated 
into 3 distinct points Q a , 1 < a < 3. Under the involution er, the leading terms in the three 
solutions ( p.4|) transform as 

(x, 1 + y x x + • • •) -> (-x, (1 - y x x + • • -) _1 ) = (-x, 1 + yis + • ■ •) (2.6) 

Since the three solutions yi of the equation fl2.5|) are distinct for generic values of the moduli 
Uj, we see that each of the three points Q a above x = are fixed under the involution er. 

• For generic values of the moduli Ui, there are also three distinct branches of y(x) near 
x = oo. A first branch y(x) = 0(x N+2 ) with a pole of order N + 2 can be readily found 

= X N+2 (u N + MAr_iX _1 + U N - 2 X~ 2 + •••)■ ( 2 -^) 

(The first three coefficients in y(x) turn out to be exactly the first three coefficients un,un-i 
and in the polynomial P/v(x) of (|2.2|) .) We denote by Pi the corresponding point above 
x = oo. In view of the involution a, a second branch y(x) = 0(x~( N+2 ^) with a zero of order 
N + 2 exists which is the image of the first branch under er 

y(x) = (-x)-( N ^±- ( 1 + ^or 1 + u n-i-vnu N -2 x _ 2 + \ 

U N \ U N U 2 N J 

The corresponding point above x = oo is denoted P3. Finally, the involution er implies that 
the third branch y(x) is regular and fixed under er 

y(x) = (-) N+2 [l + oQ] (2.9) 

Denoting the corresponding point above x = 00 by P 2 , we have 

cr : P x P 3 , o : P 2 +4 P 2 . (2.10) 

• The branching points of T over x-plane are just the zeroes on T of the function d y R(x, y) 
which are different from the singular points Q a . This function has a pole of order 2(N + 2) 
at Pi and a pole of order (N + 2) at each of the points P 2 and P3. Therefore, it has 4iV + 8 
zeros. At each of the points Q a the function d y R(x, y) has zeros of order 2. Hence 

#{Branch Points} = AN + 2. (2.11) 

Note that for generic moduli «j, neither nor 00 is a branch point, in view of our previous 
discussion. Also for generic «j, we can assume that the ramification index at all branch points 
is 2. Thus the total branching number is just the number of branch points. Since the number 
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of sheets is 3, the Riemann-Hurwitz formula can be written as g(T) = -3 + |(4AT + 2) + 1 
in this case. Thus the genus g(T) of the curve V is 

g{T) = 2N-l. (2.12) 

• For generic moduli Uj, the involution o : T — > T has exactly four fixed points, namely 
the three points Q a above x = and the point Pi above x = oo. That implies that the 
factor-curve T/a has genus 

<7(r/<7) = iV-l. (2.13) 

The involution a induces an involution of the Jacobian variety J(T) of T. The odd part 
J Pr (T) of J(T) is the Prym variety and the even part is isogenic to the Jacobian J(T/a) of 
the factor-curve T/a. The dimension of the space of divisors [D] which are even under a is 
equal to dim J (T/a) = N — 1. 



3 The Spin Models 

We introduce two systems with the same family of spectral curves fl2.1|) . One system has 
non-trivial dynamics along the even while the other system has non-trivial dynamics along 
the odd (Prym) directions of the Jacobian. The system corresponding to the SU(N) Yang- 
Mills theory with a hypermultiplet in the anti-symmetric representation is the even system. 
We sketch here the outline of the construction of both models, leaving the full discussion to 
sections §4-5. 

Both models are periodic spin chain models, with a 3-dimensional complex vector at 
each site. We view three-dimensional vectors s as column vectors, with components s a , 
1 < a < 3. We denote by s T the transpose of s, which is then a three-dimensional row 
vector, with components s a . In particular, s T s is a scalar, while ss T is a 3 x 3 matrix. Since 
the odd divisor spin model is simpler, we begin with it. 



3.1 The Odd Divisor Spin Model 

The odd divisor spin model is a (N + 2)-periodic chain of complex three-dimensional vectors 
s n — $N+n+2, s n = (s n>a ), a = 1,2,3, subject to the constraint 

3 

s n s n ^ s n s n ^ a 0, (3.1) 

a=l 

and the following equations of motion 

s n+l S n-1 cy\ 

S-n+lSn S n _iS n 

The constraint ( |3.1[ ) and the equations of motion are invariant under transformation of the 
spin chain by a matrix V satisfying the condition V T V = I 

s n -> Vs n (3.3) 
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The odd divisor spin model is integrable in the sense that the equations of motion are 
equivalent to a Lax pair. To see this, we define the 3x3 matrices L n (x) and M n [x) by 



L n (x) = l + xs n sl (3.4) 

M n (x) = X^- (Sn-l«n + s nSn-l) ( 3 - 5 ) 

S n S n -i 

Then the compatibility condition for the system of equations 

ip n+ i = L n (x)if) n (3.6) 

i) n = M n (x)^ n (3.7) 

is given by 

L n (x) = M n+1 (x)L n (x) - L n (x)M n (x) (3.8) 

A direct calculation shows that for L n (x) and M n (x) defined as in ( |3.5|) , this equation is 
equivalent to the equations of motion ( |3.2| ) for the spin model. Define now the monodromy 
matrix L(x) by 

AT+l 

L(x) = L N+1 (x) ■ ■ ■ L {x) = J] L n (x) (3.9) 

n=0 

where the ordering in the product on the right hand side starts by convention with the lowest 
indices on the right. Then L(x) and M(x) = Mq(x) form themselves a Lax pair 

L(x) = [M(x),L(x)} (3.10) 

This is easily verified using fl3.8|), since 



N+l N+l k-1 

= E II L n (x) L k x l[ L n (x) (3.11) 

fc=0 n=k+l n=0 

N+l N+l fe-1 

= E II L n (x)(M k+1 L k - L k M k ) Y[L n {x) (3.12) 

fe=0 n=k+l n=0 

N+l N+l k N+l N+l fe-1 

= E II L n (x) M k+1 H L n (x) - H L n (x) M k Y[ L n (x) (3.13) 

k=0 n=k+l n=0 k=0 n=k n=0 

= M N+2 L{x) -L(x)M {x). (3.14) 

In particular, the characteristic equation of L(x) is time-independent and defines a time- 
independent spectral curve 

r = {(x, y); = R(x, y) = det (yl - L(x))}. (3.15) 

We assert that these spectral curves are Landsteiner-Lopez curves ( |2.1| ). In fact, it follows 
immediately from the expression ( |3.5|) that detL n (x) = 1, L n (x) = L n (x) T , and L n (a;) _1 = 
L(—x). Thus 

detL(x) = 1, L(x)- 1 = L(-x). (3.16) 
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These two equations imply that det(yl — L(x)) is of the form fl2.1|) for some polynomial f(x). 
To obtain the expression ( |2.2|) for f(x), it suffices to observe that 



N+l 

f(x) =TiL(x) =Ti(l + xY / s n sl) + 0{x 2 ) = 3 + 0(x 2 ). (3.17) 

n=0 

Define the moduli Ui of the curve R(x, y) = as in ( |2.1D by f(x) = 3 + x 2 J2iLo UiX 1 . Then 
the correspondence between the dynamical variables s n , < n < N + 1, and the moduli Ui 
is given by 

Ui = ^ S ni S n2 S n2 S na • • • S n ._ l S ni (3.18) 

h 

where the summation runs over the set of all ordered z-th multi- indices n\ < n% < ■ ■ ■ < ni. 

To obtain the phase space of the model, we consider the space of all (N + 2)-periodic 
spin chains s n , subject to the constraint Q3.1J) , and modulo the equivalence 

Sn ~ V s n , where 

V is a matrix satisfying V T V = I. The dimension of this space is 

dim{s n }/{s n ~ Vs n } = 2N + 1. (3.19) 

Indeed, the (iV + 2)-periodic spin chains s n have 3(A^ + 2) degrees of freedom. The constraint 
( |3.ip removes N + 2 degrees of freedom, and the equivalence s n ~ Vs n removes 3 others, 
since the dimension of the matrices V with V T V is 3. A 2iV-dimensional symplectic manifold 
£odd - g Q^taingQi \yy setting 

£ odd = { Sn - u N = constant }/{s n ~Vs n } (3.20) 

On the space £ odd , the system is Hamiltonian with respect to the symplectic form defined 
by the differential dA(i) = (lnx)y, with Hamiltonian 

JV+l IT \ 

(1) ~ ~u~7 ~ ^ (s T s )(s T s 7) ( ' 

The action- variables are the periods of the differential dX(\) = — (lnsc)-^ over a basis of iV 
cycles for the curve T, which are odd under the involution a . If the curve T is viewed as a 
two-sheeted cover of T/a, these N odd curves can be realized as the N cuts along which the 
sheets are to be glued. 



3.2 The Even Divisor Spin Model 

The even divisor spin model is the Hamiltonian reduction of a periodic spin chain model 
which incorporates a natural gauge invariance. 

The starting point is a (N + 2)-periodic chain of pairs of three-dimensional complex 
vectors p n = (p n ,a), In = (Qn,a), 1 < « < 3, satisfying the constraints (|1.3|). We impose the 



equations of motion (|1.6|). As noted before, the constraints and the equations of motion are 
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invariant under the gauge transformations ( |1 .7| , [1.8D . In particular, a gauge fixed version of 
the equations of motion (|1.6|) is 



_ Pn+l Pn-1 . _ Qn+1 _ gn-1 /„ ^ 

Pn+l<?n Pn-lQn Ptln+l P n Qn-l 

This version follows from the other one by the gauge transformation 

Pn -> A„(t)p n , g n -> A^^gn, X(t) = exp J n n {t')d£\ . (3.23) 
We shall see in the next section that the system (1.6) admits a Lax representation. 



A reduced system is defined as follows. We impose the additional constraints ( |1.4| ). With 
these constraints, the spectral curves of the system are the Landsteiner-Lopez curves ( |2.1| ). 
The dimension of the phase space M. of all (q n ,p n ) subjected to the previous constraints and 
divided by the gauge group G of ( |1.7|jl.8| ), is 

dirndl = dim {(q n , Pn )}/G = 2N. (3.24) 

To see this, assume that N is even (the counting for N odd is similar). Then the constraint 
( |1.4| ) reduces the number of degrees of the (N + 2)-periodic spin chain (q n ,p n ) to the number 
3(iV + 2) of a (iV + 2)-periodic spin chain. The constraint (|1.3| ) and the gauge transformation 
( |1.7| ) each eliminates ^ + 1 degrees of freedom. Now the dimension of the space of matrices 
W satisfying W go = g^W is 5. However, in the gauge transformation (|1.8| ), the matrices 
W which are diagonal have already been accounted for in the gauge transformation ( |1 . 7| ) . 
Altogether, we arrive at the count which we announced earlier. 

The phase space {(q n ,Pn)}/G itself can be reduced further, to a lower-dimensional phase 
space defined by suitable constraints on the moduli space (w , ■ ■ • , Wat)- ^ turns out that 
there are 2 possible natural further reductions, each related to its own choice of differential 
dX and corresponding Hamiltonian structure: 

• On the (2iV — 2)-dimensional phase space defined by the constraints 

A^o = {(q n ,Pn); un = 1, un-i = 0}/G (3.25) 

the system is Hamiltonian with respect to the symplectic form defined by the differential 
dX = x—. Here we have used the same notation for the space just introduced and the 
space A^o described in the Main Theorem, in anticipation of their isomorphism which will 
be established later in §4. The Hamiltonian is given by H = utv-2 or equivalently by ( |1.15| ). 

The act ion- variables are periods of dX along a basis of N — 1 cycles Aj of T which are 
even under the involution a. (Equivalently, the A t correspond to a basis of cycles for the 
factor curve T/cr.) This is the desired integrable Hamiltonian system, corresponding to 
the M = 2 supersymmetric SU(N) Yang-Mills theory with a hypermultiplet in the anti- 
symmetric representation. 

• On the (2N — 2)-dimensional phase space A^2 defined by the constraints 

A4-2 = {(q n ,Pn)'i Uq = constant, U\ = constant}/G (3.26) 
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the system is Hamiltonian with respect to the symplectic form defined by the differential 
g?A(2) = — \y- This symplectic form coincides with the natural form 

u; = £ dj£ A dq n (3.27) 



with respect to which the system (|1.6|) is manifestly Hamiltonian, with Hamiltonian 



1 N+l 

H(p } q) = In Mtv = - 111 [(PnQn~l)(Pn-lQn)\ (3.28) 
A n=0 



The act ion- variables are the periods of the differential dXp) = — over again the even cycles 
Aj of the earlier case. 



4 The Direct and Inverse Spectral Transforms 

We concentrate now on the even divisor spin model. The main goal of this section is to 
describe the map stated in the main theorem, which associates to the spin chain (q n ,p n ) a 
geometric data (r, [D]) 

(q n ,p n )^(T,[D}), (4.1) 

The curve T is obtained by showing that the dynamical system (|1.6|) for (p n , q n ) admits a Lax 
representation L(x) = [M(x), L(x)], in which case T is the spectral curve {det (yl — L{x)) = 
0}. The Lax operator L(x) also gives rise to the Bloch function, which is essentially its 
eigenvector. The divisor D is obtained by taking the divisor of poles of the Bloch function. 
A characteristic feature of the even divisor spin model is that the equivalence class of this 
divisor [D] is even under the involution a . The map ([4.1|) descends to a map from the space 
of equivalence classes of (q n ,Pn) under the gauge group G to the space of geometric data 
(r, [D]). These two spaces are of the same dimension 2N: we saw this in (|3 . 24 ) for the 



first space, while for the second, the number 2N of parameters is due to N + 1 parameters 
for the Landsteiner-Lopez curves (including U/y and w/y-i), an d N — 1 parameters for the 
even divisors [D]. It is a fundamental fact in the theory that the map ( |4.1[ ) becomes then a 
bijective correspondence of generic points 

{q n ,p n )}/G^{(T,[D])} (4.2) 

We shall refer to the construction — > described above as the direct problem. The reverse 
construction <—, which recaptures the dynamical variables (p n , q n ) from the geometric data 
(r, [D]) will be referred to as the inverse problem. As usual in the geometric theory of solitons 
|2"2| |, it will be based on the construction of a Baker- Akhiezer function. We now provide the 
details. 



4.1 The Lax Representation 

We exhibit first the Lax representation for the system ( |1.6|) . The desired formulas can be 
obtained from a slight modification of the easier odd spin model treated in §3.1. Let p n , q n 
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be (iV + 2) -per iodic, three-dimensional vectors satisfying p^q n = 0, and define matrix- valued 
functions L n (x) and M n (x) by 

L n (x) = 1 +xq n p T n , M n (x) = x - . (4.3) 

\P n 1n-l Pn-lQnJ 

Then a direct calculation shows that the matrix functions L n (x) and M n {x) satisfy the Lax 
equation 

d t L n = M n+1 L n - L n M n (4.4) 
if and only if the vectors p n and q n satisfy the equations of motion (|1.6|). 

As before, the equation (|4.4j ) is a compatibility condition for the linear system ipn+i — 
L n (x)ip n , = M n (x)ip n - To obtain the spectral curve T, we observe that the same arguments 
as in the case of the odd spin model show that the matrix M(x) = M (x) and the monodromy 
matrix L{x) defined by L(x) = Il^o 1 L n (x) form again a Lax pair 

L{x) = [M(x),L(x)], (4.5) 

Thus the spectral curve T = {(x, y); R(x, y) = det(yl — L(x)) = 0} is time-independent and 
well-defined. We have used here the same notation R(x,y) as for fl2TT|), since the equation 
det (yI — L(x)) is indeed of the Landsteiner-Lopez form. To see this, we note that det L n (x) = 
1 and L n (—x) = L n (x)~ l . Together with the constraint ( |1.4[ ), this implies 

detL(x) = l, L(-x) =g L- 1 (x)g . (4.6) 

But we also have near x = 

N 

Tr L(x) = Tr (1 + x £ q n p T n ) + 0(x 2 ) = 3 + 0(x 2 ), (4.7) 

n=0 

so that det(yl — L(x)) is of the form ( |2.1| ). 

We observe that the expression R(x,y) = det (yl — L(x)) is invariant with respect to 
the gauge transformations (|1.7|) and ( |1.8|) . Therefore, if we write R(x, y) in the Landsteiner- 
Lopez form ( |2.1| ) with moduli Ui, the moduli it, are well-defined functions on the factor-space 
M.. In analogy with the odd spin case, Ui can be written in terms of the dynamical variables 

{Pn, Qn) as _^ 

u k = E(pJ?M)(P^?i3) ■ ■ • (pjfti) (4-8) 

h 

Here the summation is again over sets Ik of multi-indices J = {i\ < i 2 < ■ ■ ■ < ih)- 
4.2 General Properties of Bloch Functions 

The points Q = (x, y) of the spectral curve Y = {(x, y); det(yl — L(x)) = 0} parametrize the 
Bloch functions {ip n (Q)}o<n<N+i of the spin model. We begin by recalling the definition of 
Bloch functions, and by describing their main properties in the case of our model. 
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• We fix a generic choice of moduli parameters Ui. Then the matrix L(x) has 3 distinct 
eigenvalues y, except possibly at a finite number of points x. Let Q = (x,y). The Bloch 
solution ip n (Q) for the spin model {L n (x)}o< n <N+i is the function if) n (Q) with the following 
properties 

*n+l(Q) = L n (x)* n {Q), *JV+n +2 (<5) = J/*n(Q)- (4-9) 

These equations determine ip n {Q) only up to a multiplicative constant. To normalize tp n (Q), 
we observe that for generic moduli parameters itj, there are only finitely many points Q where 
the eigenvector V'o(Q) of the matrix L(x) satisfies the linear constraint S„ = i'0o(Q) = 0- 
Outside of these points, we can fix ip n {Q) by the following normalization condition 

X>? = 1- ( 4 - 10 ) 

a=l 

The Bloch function ip n {Q) is then determined on the spectral curve T outside of a finite 
number of points, and hence uniquely on T. Furthermore, the components of ip n {Q) are 
meromorphic functions on V. This follows from the constraint ( |4.10| ) and the equation 



L(x)i/j (Q) = yip (Q). They imply that ipo(Q) is a rational expression in y and in the entries 
of the matrix (L a p{x) — L a3 (x)) i< a <3 , in view of Cramer's rule for solving inhomogeneous 

systems of linear equations. Since x, y and L a p(x) are all meromorphic functions on T, our 
assertion follows. 

• The exceptional points excluded in the preceding construction of Bloch functions are 
the points where L(x) has multiple eigenvalues, and the points where the eigenvector ipo{Q) 
lies in the linear subspace of equation J2a=i i>o(Q) — 0. By restricting ourselves to generic 
values of the moduli Ui, we can make the convenient assumption that these two sets of points 
are disjoint. In this case, it is evident that at points where Y^i^oiQ) — 0, the function 
ipo(Q) develops a pole. 

Consider now a point Xq ^ where the matrix L(x) has a multiple eigenvalue. Let 
[X — Xq ) 1/b be the local holomorphic coordinate centered at the points Q lying above x , 
where the branching index b can be either 1 or 2. (We can exclude the possibility b = 3 by a 
genericity assumption on the moduli m.) The holomorphic function y on the surface T can 
be expanded as 

y = yo + m(x - x ) 1/b + o(x - x ), (4.11) 

where e b = 1 is a root of unity. If b = 1, it follows that 

d x R(x , y ) = d y R(x , y ) = 0, (4.12) 

which means that the curve is singular at (xq, yo). By a genericity assumption on the moduli 
Ui, the only singular point on T is at Xo = 0, and this possibility has been excluded. Thus 
b = 2, and the curve T has a branch point at xq if and only if L(xq) has multiple eigenvalues. 
The matrix L(xq) can now be shown to be a Jordan cell, i.e., L(xq) is of the form 



L(x ) = A 2 (4.13) 
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in a suitable basis, for some \x 7^ and Ai = A 2 7^ A3. In fact, L(x ) has only one double 
eigenvalue by genericity assumptions on Ui. The three branches of the function y consist 
then of one branch which is of the form A3 + y\{x — Xq) + • ■ ■ and is holomorphic in the 
variable The other two branches are of the form 

y = X 1 ±y 1 (x-x ) 1 / 2 + ---. (4.14) 

We must have y\ 7^ 0, for otherwise y = 0(x — Xq), and the same argument which ruled out 
the branching index 6 = 1 would imply that T is singular at xq. Now for x near but distinct 
from xq, the Bloch function ipo{Q) also has 3 distinct branches. Let ip± be the branches 
corresponding to the eigenvalues in (|4.14j ), and expand them as 

$ ± = </4 0) + {x - x ) 1/2 ij¥ + 0(x - x ) (4.15) 
Up to 0( X — Xq ), the eigenvector condition can be expressed as 

L(x)0/4 O) + (x- X ) 1/2 V4 1} ) = (Ai ± yx(x - x ) 1/2 )(^i 0) + (x- xo) 1 ^) (4.16) 
This is equivalent to 

L(* )V4 0) = A^f, (L(x ) - X 1 )^2 ) = ±yi^± ] - (4.17) 

Clearly, this equation admits no solution if L(x ) is diagonal. Thus L(x ) is of the form ( |4. 1 3| ) 
with /i ^ 0. We can now identify the coefficients ijj±^ and ip± in the Puiseux expansion 
( [4.15| ). The eigenspace of L(xq) corresponding to the eigenvalue Ai is one-dimensional and 
generated by a single vector 0i, which we can take to satisfy the normalization condition 
( |4.10| ). Evidently, = <pi. Let <p2 be the second basis vector in the basis with respect 
to which L(x ) takes the Jordan form ( |4.13| ), i.e., L(x )4>2 = y\4>2 + ^4>i- Then the second 
equation above is solved by 

V4 1} = ±(-0 2 + Wi), (4.18) 

II 

where the constant v is chosen so that J2a=i V'i = 0- 

• Outside a finite number of points x, the matrix L(x) has 3 distinct eigenvalues y(a) 
and three distinct eigenf unctions ^0(0), 1 < a < 3, normalized uniquely by the condition 
( |4.10| ). The function 

det 2 {^ (l) ^o(2) ^o(3)} (4.19) 

is independent of the ordering of both ipo{o) and the corresponding eigenvalues y(a). By the 
preceding observations, it can be expressed as a rational function of x and y(a), which is 
also symmetric under permutations of y(a). Thus it is actually an unambiguous and rational 
function of x. We observe that the function det 2 {ijj (l) ipo(2) "00(3)} vanishes at exactly 
those values of x which are branch points for the spectral curve det(yl — L(x)) = 0. Indeed, 
we saw earlier that the branch points Xq are exactly the points where L(xq) has multiple 
eigenvalues. Outside points xo where L(xo) has multiple eigenvalues, the determinant ( |4.19|) 
is readily seen to be 7^ (it may be infinite, because of the normalization ( |4.1(J| )). Conversely, 



assume that Xq is a branch point. Then our preceding discussion shows that for x near xq 
det 2 {^ (l) <M2) ^o(3)}(x) = (x - x o )det 2 {0 1 2 </> (3)} + 0(x - x f' 2 (4.20) 
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This shows that det 2 {^ (l) Vo(2) V>o(3)}(x ) = linw det 2 {V> Q (l) ^b(2) Vo(3)}(x) = 0, 
establishing the observation. Furthermore, since the vectors 0i, <p2, and ^o(3) are linearly 
independent by construction, we obtain the important fact that the order of vanishing of the 
square of the determinant in ( |4.19| ) at a branch point is exactly 1. (More generally, for an 



arbitrary branching index b, the order of vanishing of the square of the determinant is equal 
to b — 1, although we do not need this more general version here, thanks to our genericity 
assumption on the moduli 

• We can now determine the number of poles of the Bloch function ipo{Q) outside of the 
points P a above x = oo. Clearly, this number is half of the number of poles of the expression 
( |4.19| ) outside of x = oo. Now at x = oo, we saw that the operator L(x) has 3 eigenvalues, 



so that ( 4.19|) does not vanish there. Furthermore, we shall show later that ipo{Q) is finite at 



all three points above x = oo. Thus (|4.19|) has neither a zero nor a pole at x = oo. In view 
of the preceding discussion, the number of zeroes of (|4.19 ) is equal to the number of branch 



points of T. We showed earlier, using the Riemann-Hurwitz formula, that the number of 
branch points of F is AN + 2. It follows that the number of poles, and hence of zeroes of 
i[>o(Q) on T is 2JV + 1. 

• The poles of ip n {Q) outside of the points P a lying above x = oo are independent of n. 
To see this, we let Sq(x) be the 3x3 identity matrix J, and set S n (x) = L n -i(x)S n -i(x) = 
Ilo<fc<n-i Lk(x). Then ip n {Q) can be expressed as 

V>n(Q) = £ n -l(x)Vn-l(<2) = LI L k (x)MQ) = Sn(x)Mx) (4-21) 

0<fc<n-l 

This shows that the poles of ip n (Q) outside of P a can only occur at the poles of ipo{Q). For 
generic values of the moduli Ui, we can assume that all the poles of i> n (Q), < n < N + 1, 
are exactly of the same order 1 when they occur outside of the points P a above x = oo. 

• Let D = {zx, ■ • • , Z2N+1} be the divisor of poles of the Bloch function ip n {Q). Then a 
fundamental property of the even divisor spin chain model is the invariance of the equivalence 
divisor class [D] of D under the involution a 

[D] = [D a ] (4.22) 

In other words, there exists a meromorphic function on V with poles at z n and z*. This is a 
consequence of how L(x) transforms under the involution x — > —x, y — > y^ 1 

L(-x)=g L(x)- 1 g . (4.23) 

This transformation rule implies that goipo{Q) is a Bloch function at (— x, y~ x ) if ipo(Q) is a 
Bloch function at (x,y). Thus ^oV'o(Q) must coincide with ^(Q* 7 ) up to normalization 

9oMQ) = f(Q)MQ a ) (4-24) 

Since both ipo{Q) and 4 ! o{Q a ) are meromorphic functions, the function f(Q) is meromorphic. 
This proves ( |4.22| ). 

We summarize the discussion in the following lemma: 
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Lemma 4.1 The vector-function ip n {Q) is a meromorphic vector-function on T. Outside 
the punctures P a ( which are the points of V situated over x = oo) it has g + 2 = 2N + 1 
poles {zi, . . . , z 2 n+i}, which are n-independent. The divisor class [D] of D is invariant with 
respect to the involution a, i.e. there exists a function f(Q) on F with poles at Zj and zeroes 
atzj. 



4.3 The Direct Problem 



In the previous discussion, we made use only of the fact that the curve R(x, y) — is the 
spectral curve of a matrix L{x) which satisfies the involution condition L(—x) = L(x) . 
In particular, the discussion applies for generic values of the moduli Ui parametrizing the 
curves. 

We consider now the direct problem for the system ( |1.6| ), where the matrix L(x) arises 
more specifically in terms of the dynamical variables (q n ,p n ) as L(x) = Y\n=o L n (x) = 
Y\n=o(^ + x QnPn)- The discussion in the previous section has provided a precise description 



of the right hand side of the map (4.1). It is also evident that the map descends to the 
equivalence classes of (q n ,p n ) under the gauge group G. 

It is convenient to exploit the gauge transformation (|1.8| ) to normalize the Bloch functions 
at x = 0. First, we observe that L n (0) = I for all n, so that ip n (Qa) is independent of n. 
Furthermore, the Lax operator L(x) can be written near x = as 

L(x) = I + xT + 0(x 2 ) (4.25) 

where the matrix T is given by 

N+l 

T = £ q n pl (4.26) 

n=0 

In particular, T satisfies the condition 

T = g Tg . (4.27) 

in view of the constraint p n = g P-n-i, q n = goq-n-i- Next, recall from our discussion of the 
Landsteiner-Lopez curve in §2 that T has 3 distinct eigenvalues yi(Q a ), and that y can be 
expanded as y — l+yi(Q a )x+0(x 2 ) near Q a . Expanding ipo(Q) near Q a as ipo(Q) = ipo(Qa) + 
0(x), and using the preceding expansion for L(x), the condition L(x)ipo(Q) = yipo(Q) for 
Bloch functions can be rewritten as 

(/ + xT)(MQa) + x^' (Q a )) = (1 + y a x)(MQa) + xipoiQa)) + 0(x 2 ) (4.28) 
This implies 

TMQa) = VaMQa) (4.29) 

i.e., ipo(Q a ) are precisely the three eigenvectors of T, corresponding to the eigenvalues y a . If 
we let \J/ (0) be the 3x3 matrix whose columns are the vectors ipoiQa), then the transfor- 
mation law (|4.27| ) implies that \J/ (0) satisfies the condition 

*o(0) = g ^o(0)g . (4.30) 
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Now the transformation ( |1.8| ) on (q n: p n ) does not change the curve V and the divisor D, but 
changes the matrix \I/o(0) into W^q(0). But \I/o(0) commutes with the matrix g , and hence 
so does its inverse. This means that the inverse qualifies as one of the gauge transformations 
W allowed in ( |1.8|) . Under such a gauge transformation W, the Bloch function \l/o(0) gets 
transformed to the identity 

*o(0) = /• (4.31) 
Henceforth we can assume then this normalization, and p n , q n satisfies the condition 

N+l 

Ti = E Qn,y n = y?C (4.32) 

n=0 



Our main task is to establish that the map ( |4.2| ) is generically locally invertible. This is 
the goal of the next section on the inverse spectral problem, but in order to motivate the 
constructions given there, we identify here the basic behavior of the Bloch function if) n (x,y) 
near the points P a above x = oo. For (x,y) near P a , set 

oo 

V„(z, y) = x p - £ ^ n ±{P^~ k (4-33) 

k=0 

Here p na is the order of the pole (or zero when p na < 0) of ip n (x,y) near P a , which may 
vary with both n and a. The following lemma identifies the coefficients ip n ,k{Pa) up to 
normalization: 

Lemma 4.2 • In the neighborhood of the puncture P\ (where y = 0(x N+2 )), the vector- 
function ifj n has a pole of order n and the leading coefficient ip n ,o{Pi) of its expansion is equal 
to 

1pn,o( P l) = a nq n -l, (4.34) 

where the scalar a n satisfy the recurrence relation 

a n+1 = (plq n -i)a n . (4.35) 
The next coefficient ^ n ,i(Pi) satisfies 

lpn+1,1 = ll>n,0 + QniPn^n,!) (4.36) 

• In the neighborhood of the puncture P3 (where y = 0(x~ N ~ 2 )) the vector-function ip n has 
a zero of order n and the leading coefficient ipn,o{P^) of its expansion is equal to 

i>n,o{ P 3) = Pnq n , (4.37) 
where the scalar (3 n satisfies the recurrence relation 

Pn + 1 = -T^—^n- (4-38) 
[PnQn+l) 

• In the neighborhood of the puncture P2 (where y = 1) the vector-function ip n is regular and 
its evaluation ^,0(^2) at P2 is orthogonal to both p n and p n -i, i-£-, 

p^n,Q(P 2 ) = pLiVVo^) = 0. (4.39) 
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Proof. First, we show that for generic moduli, the Bloch function ipo(x, y) is regular near each 
P a . Observe that ip^ +2 {x,y) = L(x)ip (x,y) = yip (x,y). Now the relation ip n+ i = L n (x)^ n 
can be inverted to produce 

ip n (x, y) = L n (x)"Vn+i(a;, y) = (1 - xq n pl)ip n+1 (x, y) (4.40) 

Applying this relation N + 2 times, we may write 

ipo(x, y) = 2/(1 - xqoPo ) • • • (1 - xq N+ ip T N+l )i) {x, y) (4.41) 

Consider first the neighborhood of the point P3, where y is of order x~^ N+2 \ If ?po(x,y) 
admits the expansion (|4.33| ) near P 3 with V'o^-F's) 7^ 0, then we must have 

^0,0 = (-f +2 qo(Po qi) • • ■ (pW+iKPn+iM (4.42) 

This shows that ^0,0 (-P3) is proportional to the vector q , say ^0,0 = Pqo- Now recall that 
the Bloch function ip (x,y) satisfies the normalization condition ( |4.10| ) throughout. This 
implies that J2a=i ^0,0(^3) = if the order n (P 3 ) of the pole of i[) (x,y) at P 3 is positive. 
For generic values of the moduli of the curve T, we may assume that Z)a=i ^oa 7^ 0. It 
follows that Po = and hence ^0,0 (-P3) — 0, which contradicts the definition of ^0,0(^3)- 
This shows that no(P^) = 0, and the Bloch function ipo(x, y) is regular at P3. The arguments 
near Pi is similar and even more direct, just using the equation yip (x,y) = ipN+2(x,y) = 
Tln=o(l + zqnPn)^>o(x, y) ■ It shows, incidentally, that the leading coefficient Vo,o(-Pi) is 
proportional to gjv+i- At P 2 , the regularity of ^0(2;, y) follows from the regularity of ipo(x, y) 
at the other two points P\ and P3, and from the fact that for generic moduli, the determinant 
( |4.19| ) is regular. 

It is now easy to see that the functions ip n (x, y) have the zeroes and poles spelled out in 
Lemma 4.2. The recurrence relations stated there can also be read off the defining relations 
ipn+i = L n (x)ip n (x) . For example, near P 1; we find 

x n+1 Un+1,0 + ^n+1,1 + •••)= x n {l + xq n p T n ) U nfi + -^,1 + • • ■) (4.43) 
This implies 

= qnifn^nfi) (4.44) 

lpn+1,1 = VVO + qn(Pn^n,l) (4.45) 

The relations (|4.35| , |4.36|) follow. Near P 2 , we write instead 

X~ n U n ,0 + -^n,l + •••)= 2T n-1 (l - Xq n pl) U n+ i, + ~V>n+l,l + (4-46) 

This implies 

^nfl = -g n (Pn^n+l,o) (4.47) 

ip n>1 = tpn+i,o ~ g„(^Vn+i,i) (4-48) 
which gives ( |4.37| , 4.38 ). Finally near P 2 , we get 



i> n +i,o + -VVt-1,1 + 
x 



(l + xq n pl) (^ n ,0 + ^n,l + • • •) (4.49) 



This implies that p^Vyo = 0. Furthermore, VWi,o — V^o + qnip^^n,!) ■ Multiplying on the 
left by we conclude that p^ipn+ifi = 0. This establishes ( |4.39| ), and Lemma 4.2 is proved. 
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4.4 The Inverse Spectral Problem 



It is now a standard procedure in the geometric theory of soliton equations to solve the 
inverse problem using the concept of the Baker-Akhiezer function originally proposed in 



22]. The main properties of the Baker-Akhiezer function in our model are the following. 



• Let T be a Landsteiner-Lopez curve defined by equation ( j2.1|) . Then for a divisor D 
of degree g + 2 = 2N + 1 in general position, there exists a unique vector-function <p n (t, Q) 
such that: 

(a) 4> n (t,Q) is meromorphic on T outside the punctures Pi,P 3 . It has at most simple 
poles at the points Zi of the divisor D (if all of them are distinct); 

(b) In the neighborhood of the punctures Pi and P3, it has respectively the form 



4> n = x n e xt X>»,*(Pi)aT* , Q -»• Pi, (4-50) 

\k=0 ) 

n = x- n e~ xt ff; 0„, fc (P 3 )x- fc ) , Q ^ P 3 , (4.51) 

(c) At the points Q a , 4> n (Q) is regular, and 4> n (Qa) is equal to 

(t>n,a{t, Qp) = 5 a>p . (4.52) 

The arguments establishing the existence of the Baker-Akhiezer function <p n are well- 
known, so we shall be brief. First, we recall that as shown in [T2IJ for any algebraic curve with 
two punctures, any fixed local coordinate in the respective neighborhoods of the punctures, 
and for any divisor D of degree g there exists a unique (up to a constant factor) function 
with the analytic properties stated above. Now let (Pi,P 3 ) be the punctures, and let x' 1 
be the local coordinate near either one of the punctures. We can easily show that if D 
has degree g + 2, the dimension of the space of such functions is equal to 3. We form the 
3-dimensional vector whose components are just the three independent functions from this 
space. This 3-dimensional vector is unique up to multiplication by a constant matrix. We 



fix this matrix by the normalization condition ( 4.52 ). This establishes our claim. 

The function <p n (t,Q) can be written explicitly in terms of the Riemann ^-function as- 
sociated with r. The ^-function is an entire function of g = 2N — 1 complex variables 
z = (zi, . . . , z g ), and is defined by its Fourier expansion 

9( Zl , ...,z g )= J2 mezg e 2m<m ^ >+m<Tm > m> , 

where r = Tjj is the period matrix of T. The ^-function has the following monodromy 
properties with respect to the lattice Z 9 + rZ 9 

0(z + /) = 8(z), 9(z + rl) = exp[-i7r <rl,l > -2m < I, z >] 9{z) 

where / is an integer vector, / e Z 9 . The complex torus J(T) = C 9 /Z 9 + rZ 9 is the Jacobian 
variety of the curve T. The Abel transform 

rQ 



T3Q^ A k (Q) = f doo k 
JQo 



19 



imbeds the curve T into its Jacobian variety. Here dwk is a basis of g holomorphic differentials, 
normalized as dual to the A-cycles of a symplectic homology basis for T. 

According to the Riemann-Roch theorem, for each divisor D — z\ + . . . + z g+ 2 in the 
general position, there exists a unique meromorphic function r a {Q) with r a {Qp) = 5 a p and 
D as the divisor of its poles. It can be written explicitly as (see details in 



where 



UQ) , fn s_ 0(A ,n\ ■ 7 ^IW»W)+ f i>) 



UQ) = "FT^T, UQ) = 0{A{Q) + Z, 



UQaY Ja ^' v ^' 'Y[ l m=l 0(A(Q) + S n 



".])■ 



9-1 9-1 

F p = - K, - A(Q P ) - A&j), S m = -JC- A(z g _ 1+m ) - £ A( 

3=1 3=1 
S+2 3 

Z a = Z - A(R a ), Z = -JC-J2 A(zj) + J2 A(Q a ), 

j=l a=l 

where K, is the vector of Riemann constants. 

Let dQo and dfl\ be the unique normalized meromorphic differentials on T, which are 
holomorphic outside Pi and P3, and with the property that dflo has simple poles at the 
punctures with residues =pl, dQi is regular at P 3 , and has the form dQi = dx(l + 0(x~ 1 )) 
at Pi. The normalization means that the differentials have zero periods around A-cycles 

f d£lo = f dCli = 0. 

J A J A 

We observe that the differential dQl(Q) = dQ(Q a ) has a pole only at P3, and is there of the 
form — dx(l + 0(x~ l ). 

Let V and U be the vectors whose components are the P-periods of the differentials dfl 
and dQi respectively 

v = ^-i dQi, u = ^-l dn . 

2m Jb 2tti Jb 

The Baker- Akhiezer function (p n (t,Q) is given by 

Q) = UQ) e{A{Q) + Za) e(tu+ + nv + z ) exp [L ndn ° + tdn ) (4 - 53) 

where dtt + = dQ x + dQf and U + = U + U a . 

• The Baker-Akhiezer function n is a Bloch function, in the sense that 

(PN+2+n (t,Q) = y<f> n (t,Q). (4.54) 

This is just a consequence of the fact that both sides of the equation satisfy the criteria for 
the Baker-Akhiezer function, and that the Baker-Akhiezer function is unique. Similarly, the 
uniqueness of the Baker-Akhiezer function implies that, if the divisor D is equivalent to D a , 
then the function (p n satisfies 

Mt,Q)=9o<f>-n(t,Q a )f(Q), (4-55) 
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where f(Q) is a function with poles at 7 S and zeros at 7J. Without loss of generality, we 
may assume that f(Qi) = /(Q3) = — /(Q2) = 1- 

• Let 4> n {t,Q) be the Baker-Akhiezer function corresponding to T and the divisor D 
of degree 2g + 2. Let p„,(t) be a vector orthogonal to <p n ,o(P3,t) (the leading term in the 
expansion ( f4.50| )), and to (f> n (t, P2), i.e. 

^0 n , O (P 3 )=Pn0n(t,P 2 )=O, (4.56) 

and q n be the vector 

P»0n-l,o("l) 

The vector functions p n , q n are then (N + 2)-periodic and mutually orthogonal, and they 
satisfy the contraint ([4.31|) . As can be expected from the gauge invariance ( |1.7| ) in the 
direct problem, the functions p n {t) and q n {t) which we obtain this way are defined only 
up to a multiplier /z n (t). However, the operators L n and M n (x) are uniquely defined by 
the expression ( f4.3| ). Furthermore, again by uniqueness of the Baker-Akhiezer function, the 
Baker-Akhiezer function <f) n (Q) satisfies 

Vvh-i(*> Q) = L n {x)Mt, Q), (dt - M n {x))tfj n {t, Q) = 0. (4.58) 

Thus the vector function (q n (t) , p n (t) is a solution of the dynamical system (|1.6| ). If the 
equivalence class of the divisor D is invariant with respect to a, then (p n , q n ) satisfies in 
addition the relation (|1.4p. 

• The Baker-Akhiezer function <p n (t, Q) satisfies the same defining Bloch property ( |4.54|) 
as the Bloch function ip n (Q), except for the different normalizations, which is ( |4.52| ) in the 
case of <fi n (t, Q) and ( [4.10|) in the case of ip n (Q)- It follows that 

3 

Mt,Q) = r'\t,Q)(P n (t,Q), r(t,Q) = Y,<f>o,a(Q) ( 4 -59) 

0=1 

is a Bloch solution of ( |4.9| ) normalized by the condition ( |4.10[ ). This leads to the following 
description of the dynamical system (|1.6|) . 

Let p n (t),q n (t) be vector functions (subject to constraints ( |1.3|J1.4| , P:.31| ) ) which satisfy 
the equations (|1.6| ). Then the t-dependence of the divisor D under the map ([4.1|) 



2Af+l 

(p n (t), q n (t)) {r, D(t) = J2 ( 4 -60) 

i=i 

coincides with the dynamics of the zeroes of the function r(t,Q) given by ( |4.59| ). The 
dynamics of the Bloch eigenfunction of ( f4.9|) (i.e. normalized by ( [4.10|) ) are described by 

(dt-M n (t,x))4> n (t,Q) = ii(t,Q)il> n (t,Q), fi = -d t \nr(t,Q) (4.61) 

We observe that the linearization of the equations of motion on the Jacobian of the curve is 
a direct corollary of the linear dependence on t of the exponential factor in the expansion of 
ip n (t, Q) near the punctures. 
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As we saw earlier, the normalization (|4.31) can be achieved by the action (|1.8|) of a 



subgroup of matrices W which commutes with g . In order to get A4, we have to consider 
in addition the action of diagonal matrices W. The basic observation is the following. 

Let the Baker- Akhiezer functions if) n (t, Q) and ip'(t, Q) corresponds to equivalent divisors 
D and D', respectively. Then 

H(t,Q) = WMt,Q)KQ), (4-62) 

where h(Q) is a function with poles at D and zeros at D', and W is a diagonal matrix 
W = h~ l {Q a )8 a fi. To establish ( }4.62|) , it suffices to check that both sides of the equation 



have the same analytical properties. The equation ( [4.62| ) implies that the vectors p n , q. 



defined by equivalent divisors are related by a transformation (|1.8|) with a diagonal matrix 
W. Altogether, we have established the following part of the Main Theorem of §1: 

Theorem 1. The map identifies the reduced phase space M. with a bundle over the 

space of algebraic curves Y defined by ( \2. 1\ ) with /at(x) of the form ( \2.2\ ). At generic data, 
the map has bijective differential. The fiber of the bundle is the Jacobian J(T Q ) of the factor- 
curve r = r/cr. 

M = {T,[D]e J(r )} (4.63) 



5 Hamiltonian theory and Seiberg-Witten differential: 
The Even Divisor Model 

We come now to the crucial issue of how to determine the symplectic forms with respect 
to which the system ( |1.6| ) is Hamiltonian. For this, we rely on the Hamiltonian approach 
proposed in || and for general soliton equations expressible in terms of Lax or Zakharov- 
Shabat equations. This approach was effective in the study of gauge theories with matter in 



the fundamental representation. Further applications were given in |T0[ and |24| • We review 
its main features. 



5.1 The Symplectic Forms in terms of the Lax Operator 

In order to find the Hamiltonian structure of the equations starting with the Lax operator, 
we need to identify a two-form on the phase space M. of vectors (<?„,£>„), written in term of 
the Lax operator L(x). Candidates for such two-forms are 

1 3 dx 
"V) = 9 E Res Pa < % +1 (Q)8L n (x) A 6* n (Q) > — (5.1) 

Z a=l X 

The various expressions in this equation are defined as follows. The notation < f n > stands 
for a sum over one period of the periodic function f n : 

N+l 

< fn >= E U (5-2) 
n=0 
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The expression tp^(Q) is the dual Baker- Akhiezer function, which is the row- vector solution 
of the equation 

r n+1 (Q)L n ( Z ) = r n (Qi r N+2 {Q) = (5.3) 

normalized by the condition 

r (Q)MQ) = i- (5-4) 

Note that (jO|) and Q imply that V£+iVVh-i = ip* +1 (L n {x)i^ n ) = (ip^ +1 L n (x))?p n = ^ip n 
does not depend on n. We would also like to emphasize that, unlike the Bloch function 
ipniQ) which does not have n- independent zeroes, the normalization ( |5.4| ) allows the dual 
Bloch function ip^iQ) to have such zeroes. In fact, they occur at the poles of ip n {Q). 

In (|5 . 1|) , the differential 5 denotes the exterior differential with respect to the moduli 



parameters of M.. (This is in order to distinguish S from the differential d, which is the 
exterior differential on the surface T.) Thus the external differential 5L n (z) can be viewed 
as a one-form on A4, valued in the space of operator- valued meromorphic functions on T. 
Similarly the Bloch function ip n {Q) and dual Bloch functions i[>n(Q) are functions on A4, 
valued respectively in the space of column-vector-valued and the space of row- vector- valued 
meromorphic functions on T. It follows that 8ijj n (Q) is a one-form on A4, valued in the space 
of column-vector-valued meromorphic functions on T. The expression ip^ +1 5L n (x) A 5ip n (x) 
is then a two- form on A4, valued in the space of meromorphic functions on T, and for each 
m integer, the expression 

dx 

tym) =< i>n+l SL n(x) A #„(x) > — (5.5) 

is a meromorphic 1-form on T. This justifies ( |5.1| ) as a two-form on Ai. 

In ( |5.1|) , we have allowed for a later choice of an integer m. We shall see shortly that 
holomorphicity requirements restrict to < m < 2, and that the symplectic form of the 
M = 2 SUSY with a hypermultiplet in the anti-symmetric representation is obtained by 
setting m = 0. 

Sometimes it is useful to think of the symplectic form u as 

oj (m) = ^Res x . =00 Tr < (^- l +l (x)5L n (x) A 8V n {x)) > (5.6) 

where ^ n ( x ) is a matrix with the columns i/; n (Qj(x)), Qj(x) = (x,yj) corresponding to 
different sheets of T. The matrix x I / n (x) is of course not defined globally. Note that ip^iQ) 
are the rows of the matrix ^/~ 1 (x). That implies that ^(Q) as a function on the spectral 
curve is meromorphic outside the punctures, has poles at the branching points of the spectral 
curve, and zeroes at the poles Zj of \& n (Q). These analytical properties will be crucial in the 
sequel. 



5.2 The Symplectic Forms in terms of x and y 



A remarkable property of the symplectic form defined by Q5.1]) in terms of the Lax oper- 
ator L(x) is that it can, under quite general circumstances, be rewritten in terms of the 
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meromorphic functions x and y on the spectral curve Y. More precisely, we have 

2N+1 8x 

U( m ) = - Yl 6ln v( z i) A — (*•)■ ( 5 - 7 ) 

i=l X 

The meaning of the right hand side of this formula is as follows. The spectral curve is 
equipped by definition with the meromorphic functions y(Q) and x(Q). Their evaluations 
x(zi), y(zi) at the points z^ define functions on the space A4, and the wedge product of their 
external differentials is a two- form on Ai. 



The proof of the formula (|5.7| ) is very general and does not rely on any specific form of 
L n . For the sake of completeness we present it here in detail, although it is very close to the 



proof of Lemma 5.1 in |24[ . 

Recall that the expression Qf m ) defined in ( |5.5| ) is a meromorphic differential on the 
spectral curve T. Therefore, the sum of its residues at the punctures P a is equal to the 
opposite of the sum of the other residues on T. For m < 2, the differential fi( m ) is regular 
at the points situated over x — 0, thanks to the normalization ( |4.31]) , which insures that 
8ip n {Q) = 0(x). Otherwise, it has poles at the poles of ip n {Q) and at the branch points 
Si, where we have seen that if>n+i(Q) ^ as P°l es - We analyze in turn the residues at each of 
these two types of poles. 

First, we consider the poles Zi of ip n {Q)- By genericity, these poles are all distinct and of 
first order, and we may write 

4>n = 4>n,o( Z i) _ 1 / \ H ( 5 - 8 ) 

It follows that 5ip n has a pole of second order at Zi 

<ty n = tM*t) , 5X( ?\v +--- ( 5 - 9 ) 
(x - X[Zi)Y 

In view of the fact that V'n+i ^as a simple zero at Zi and hence can be expressed as 

VC+i = C+i,o( x - x &)) + ■ ■ ■ ' ( 5 - 10 ) 

we obtain 

Sx Sx 

The key observation now is that the right hand side can be rewritten in terms of the mon- 
odromy matrix L(x). In fact, the recursive relations ( f4.9|) and ( |5.3|) imply that 

/ N+l \ /n-1 \ 

< l/j* n+1 5L n ^J n > = < lP* N+2 [] L m) $ L n f II L m) ^0 > (5-12) 

\m=n+l / \m=0 / 

N+l I JV+1 \ /n-1 \ 

= II L m \SL n il[L m U (5.13) 

n=0 \ 171=12+1 J \m=0 / 

= t/)* N+2 5Lil> Q = il> Q 5liiyi/) . (5.14) 
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In the last equality, we have used the standard formula for the variation of the eigenvalue of 
an operator, ip^SLipQ = ip^SytpQ. Altogether, we have found that 



Sx 

Res 2i ^ (m) = 5\ny(zi) A — (zj). 



X 1 ' 



(5.15) 



The second set of poles of Q( m ) is the set of branching points s« of the cover. The pole of 
■0* at Si cancels with the zero of the differential dx, dx(si) = 0, considered as a differential 
on r. The vector-function ip n is holomorphic at S{. However, 5ip n can develop a pole as we 
see below. If we take an expansion of ip n in the local coordinate (x — x(sj)) 1//2 (in general 
position when the branching point is simple) and consider its variation we get 



= 4>n,0 + 1pn,±(x ~ x{Si)) l/2 + 



8A 



n,± 



Sx(st 



(x - x(Si)) 1 /' 



+ 



(5.16) 
(5.17) 



Comparing with ^ = I^ ± — -\ , we may write 

di> n 



dx 



This shows that 5^/ n has a simple pole at Sj. Similarly, we may write 

dy 



5y = --£-5x( Si )+0(l). 



(5.18) 



(5.19) 



The identities (]5.18 ) and Q5.19 ) imply that 



Res,.f2 



(to) 



Res, 



< t/j^ +1 5L n dip n > A 



Sy dx 
x m dy 



(5.20) 



Arguing as for (|5.12 ), this can be rewritten as 



ReSo.fi 



Res, 



(ip* NJr2 5Ldij)^} A -j 



5ydx 
x m dy 



(5.21) 



Due to the antisymmetry of the wedge product, we may replace 5L in ( 5.21]) by (SL — dy). 
Then using the identities 



ip* N+2 (6L - 5y) = 5ip* N+2 (y - L) 
(y - L)dip = (dL - dy)ip , 
which result from ip* N+2 (L — y) = (L — y)ipo — 0, we obtain 



Res,„f2 = Res, 



* N+2 (dL - dy)ipo) A 



5ydx 
x m dy 



(5.22) 
(5.23) 



(5.24) 
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Now the differential dL does not contribute to the residue, since dL(si) = 0. Furthermore, 
^* N +2^o = STVoV'o = V' 1 - Tnus = ~ ?T%- Exploiting again the 

antisymmetry of the wedge product, we arrive at 



Res Si Q = Res Si (^* N+2 5ipo) A by 



dx 



x" 



(5.25) 



Recall that we have normalized the Bloch function ipo(Q) at x = by ( |4.31| ), and that 
near x = 0, the function y is of the form ( p.4|) . Thus 5ipQ = 0(x) and 5y = 0(x) near x = 0, 
and the differential form 

\ dx 

'f N+2 5^ )A5y—. (5.26) 



x' 



is holomorphic at x = for < m < 2. It is manifestly holomorphic at all the other points 
of T, except at the branching points s« and the poles Z\, • • • , z 2 n+i- Therefore 



^Res Si (ip* N+2 5if) ) ASy 



dx 



x' 



2N+1 j x 

J2 Res Zi (ip* N+2 5il)o) A 5y- 



i=l 



X' 



(5.27) 



Using again the expressions ( p. 16] , |5.18|) for ^ and 5ipo, and the fact that ip* N+2 = V Vo 
the right hand side of ( |5.27| ) can be recognized as 



^ b\ny{z>) N 6 ri: ' 

i=l 



X m (Zi) 



(5.2i 



The sum of ( |5.15|) and (|5.28| ) gives ( |5.7|) , since 

2N 

2^( m ) = -^Res 2i n (m ) - ^Res Si f2( m ). 
The identity ( p.7| ) is proved. 



(5.29) 



8=1 



5.3 Action-Angle Variables and Seiberg-Witten Differential 

The expression ( |5.7|) for the symplectic form &( m ) suggests its close relation with the following 
one-form on V ^ 

dX {m) =\ny — (5.30) 

Strictly speaking, the form d\ m ) is not a meromorphic differential in the usual sense, because 
of the multiple-valuedness of In y. However, the ambiguities in In y are fixed multiples of 
27T2, which disappear upon differentiation. Thus, the form dX^ is no different from the 
usual meromorphic differentials, as far as the construction of symplectic forms is concerned. 
Also, the form c?A( m ) and the form ^jX^ m+1 (for m/1; for m — 1, —(In x)^-) differ by 
an exact differential, and we shall not distinguish between them. From this point of view, 
the Seiberg-Witten form ( |1.1|) can be identified with the form — <iA(o). 

Our spin chain model has led so far to a 2iV-dimensional phase space A4, equipped with 
several candidate symplectic forms o;( m ), 1 < m < 2. We still have to reduce M. to a 
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(2N — 2)-dimensional phase space, and to identify the correct symplectic form. Remarkably, 
both selections are tied to a key physical requirement for the one-form which corresponds 
to the Seiberg-Witten of a Af = 2 SUSY gauge theory, namely the holomorphicity of its 
variations under moduli deformations. 

It is an important feature of Af = 2 Yang-Mills theories that the masses of the theory 
are not renormalized. Since the masses of the theory correspond to the poles of the Seiberg- 
Witten differential dX, it follows that 5dX must be holomorphic. Thus we need to examine the 
poles of 5 dX = Slny-^k, and identify the subvarieties of A4 along which SdX is holomorphic. 
There are 3 such subvarieties, corresponding to the choices of m: 

• On the variety A4 2 = A4 fl {w — c o, ui = ci}, the differential 5 dX@) = (S\ay)^§ 
has no pole at Q a , since y = 1 + 0(x 2 ) near x = 0. On the other hand, the differential % 
vanishes at x = oo, so 5\ay-^ is also holomorphic there, and odX{2) is holomorphic. 

• On the variety A4 — Add {u^ = 1, u N _ x = 0}, the differential 5 dX^ = (Slay)dx 
is automatically holomorphic at x = 0. Near oo, in view of the expansion () for y, we have 
5lny = 0(x 2 ) if we vary only the moduli within M.2- Thus 5dA(o) is holomorphic. 

• On the variety Mi = M (1 {un = 1}, the differential S dA(i) = (Slny)^- is still 
holomorphic, because S lay = 0(x). Near x = oo, the sole constraint {mtv-i = 1} suffices to 
guarantee that Slay — 0(~). Thus SdX^j is holomorphic. 

When m and hence dX< m \ is even under the involution a, action-angle variables can 
be introduced as follows. Restricted to A4( m ), SdX^ is holomorphic, and hence can be 
expressed for suitable coefficients <5cij as 

2N-1 

5dX( m ) = (8<k)dui, (5.31) 
i=i 

where dui is a basis of 2^ — 1 holomorphic one-forms on T. Since dXf m ) is even, only 
holomorphic one-forms du>i which are even can occur on the right hand side. We identify 
such forms with forms on T/a. We choose a symplectic homology basis Ai,Bi and a dual 
basis of holomorphic forms dui, 1 < % < iV — 1, for the factor curve T/a. The variables 
and ani can then be defined by 



dX( m ), a Di = f dX {m) . (5.32) 

Ai J Bi 

The interpretation of the variables cij is as action variables from the viewpoint of the spin 
model and as vacuum moduli from the viewpoint of the Af = 2 SUSY gauge theory. Evi- 
dently, their variations coincide with the 5ai of the equation (|5.31| ). 



Next, the angle variables fa, 1 < i < N — 1, are defined by 



2N+1 



D = {z U - ■■,Z 2 N+l} 1 ><Pi= / duJ i ( 5 - 33 ) 
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We claim now that, for m even, the symplectic form u)( m ) is a genuine symplectic form when 
restricted to M.r m ), and that a, and fa as defined above are action- angle coordinates for u)( m ) 

N-l 

U M = ^2 5^ A 5 fa on M( m y (5.34) 
i=i 

To see this, we evaluate the two-form 5(Y<j=i 1 Jq 5 dX) in two different ways. Substituting 
in ( |5.31|) , we find that it is equal to 

JV-l N-l 

5(^2 ddifa) = ^2 54>i A5di. (5.35) 
t=i i=i 

On the other hand, we can also write 



[2N+1 z \ (2N+1 z d \ 2N+1 g , , 

J S L SiX ) =* £ I (,h,B) £ = 5 SS A( ^^- (536) 



,ti ^ m (^) 



Comparing the two formulas, and making use of ( p.7| ), we obtain the desired equation ( |5.34j ). 

We observe that for the present even divisor spin model, the space Ai\ and the form 
(i) are not applicable. In fact, there are difficulties with both the dimension of M.\ which 
is odd, and the angle variables 0j defined by ( |5.33| ), which would vanish identically because 
the class of the divisor D is even. 

For the M = 2 SUSY Yang-Mills theory with a hypermultiplet in the antisymmetric 
representation, the spectral curves are given by A4 . The symplectic form is then um\, which 
provides an independent check of the choice of Seiberg-Witten form found by Landsteiner 
and Lopez. 



5.4 The Hamiltonian of the Flow 

We show now that the even divisor spin model is a Hamiltonian system. More precisely, 
restricted to each of the phase space JA(p) or M.{2)i the system is Hamiltonian with the 
corresponding symplectic form, with a corresponding Hamiltonian. We would like to stress 
that, once again, the arguments to these ends are quite general, and use only the expression 
for uj(m) in terms of the Lax operator. 

Lemma 5.1 Let m be either or 2. Then the equations ( \1.6] ) restricted on Ai( m ) are 
Hamiltonian with respect to the symplectic form W( m ) given by ( \5.J\ ). The Hamiltonians H( m ) 
are given by 

#( ) = wtv-2 (5.37) 

JV+l i N+l 

H {2) = \I1U N = J2 HPn<ln-l) = ~ HiPiVn-ljiPn-lVn)} (5.38) 

n=0 / n=0 
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Proof. By definition, a vector field d t on a symplectic manifold is Hamiltonian, if its con- 
traction %Q t u{X) = u}(X,d t ) with the symplectic form is an exact one-form 5H(X). The 
function H is the Hamiltonian corresponding to the vector field d t . Thus 



~ Res ^ (< ^n+l SL n1pn > ~ < ljj* n+1 L n 5ljj n >) ^ 



(5.39) 



Now under the flow ( |1.6| ), the Lax operators L n (x) flow according to the Lax equation ( |4.4| ) 
while the Bloch function ip n flow according to ( [4.61 ). Consequently, 



1'dt UJ {m) 



i Res P Q (< ip* n+1 5L n (M n + n)$ n > - < i)* n+l {M n+l L n - L n M n )bi) n >) 



(5.40) 

Since L n ip n = ^ n+1 , it follows that i)* n+1 M n+1 L n 5i) n = i)* n+1 M n+1 ^ n+1 - i(;* +1 M n+ i6L n if) n . 
Upon averaging in n, we obtain 



< <4)* n+l (M n+l L n - L n M n )54) n >= - < ipn +1 M n+l 5L n il) n > 
For all n, both SL n (x) and M n (x) vanish at x = 0. The differential form 

dx 



< ijj* (5L n M n + M n+1 5L n ) ^ n > 



x 1 



(5.41) 



(5.42) 



is thus holomorphic at x — 0, in both cases m = and m = 2. As we have seen, outside 
of x = oo, the poles of tp^+i are a ^ the branch ponits and are cancelled by the zeroes of dx 
there, while the poles of ip n are cancelled by the zeroes of VCfi- Thus the above differential 
form is holomorphic outside of x = 0. The sum of its residues at P a must be zero 



dx 



x' 



Res Pa < i)* n+l (5L n M n + M n+1 5L n ) tp n > 

a 

The expression ( |5.40|) for id t ^( m ) reduces to 

-^Res Pa (< ij)* n+l 8L n ip n > fJt(Q,t) 



0. 



dx 

x™ 



Applying the arguments leading to ( 5.12 ), we obtain 



^£res Po <f(]ny)/i(*,Q)^ . 



(5.43) 



(5.44) 



(5.45) 



As follows from ( |4.50|j4.5T| ), and ( 4.61 ) the function n(t,Q) is holomorphic at Pi, while it 
has the following expansion at the punctures Pi, P3 



H(t, Q) = -x + 0(1), Q -> Pi; n(t, Q) = x + 0(1), Q -> P 3 



(5.46) 



We consider now the cases m = 2 and m = separately. When m = 2, the form //^§ is regular 
at P 2 , an d has simple poles with opposite residues at Pi and P 3 . Since 8\a.y = Sun + 0(~) 
near P 1; it follows immediately that 



id t ^(2) = 8(lnu N ) 



(5.47) 
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When m = 0, we observe that the form (5 In y)dx is regular at x — oo. Indeed, the con- 
straints un = 1, Wjv-i = defining the phase space A4 in this case imply that 5lny = O(X) 
near all three points Pi, P 2 , and P 3 . For P 1 and P 3 , this statement is a direct consequence of 
fl2.7|) and Q2.8|) . For P 2 , this follows from the fact that the three roots y a of the Landsteiner- 
Lopez curve fl2.1| ) must satisfy IIa=i Va — 1- Returning to the residues in ( |5.45|) , we see that 
the point P 2 does not contribute. As for the points Pi and P 3 , they contribute exactly the 



coefficient un-2 in the expansions (|2.7| ) and (|2.8| ) for ?/ 



«9*^(o) = 5u N - 2 - (5.48) 

The lemma is proved. 



5.5 The Symplectic Form in Terms of (p n , q n ) 

The expression fl5.1|) for the symplectic forms co>( m ) in terms of the Lax operator also provides 
a straightforward way of writing ui m ) i n terms of the dynamical variables {q n ,p n )- Such an 
expression for the form appears complicated. But it is quite simple for the form cu( 2 ), 
and we derive it here. 

We have 5L n = x5(q n p^), and the contributions of the three points P a above x = oo can 
be evaluated as follows. 

At the point Pi, y = 0(x N+2 ), ip n = 0(x n ), ip n+1 = 0(x- {n+1) ), and thus the differential 
< Vn+i^n A 5ip n > ^§ is regular. The residue at Pi vanishes. 

At the point P 2 , ip n and VC+i are regular. Using the same notation as in (j4.33|) , we write 



Tpn = 1>n,0 + VW X 1 H 

■0n+l = 0n+l,O + "0n+l,l x + " " " 

In analogy with (ref), from the equation 

C+l = CM^r 1 = Ki 1 - xq n pl), 

it follows that 

C,0& = C+l,0?n = 

The residue at P 2 is then readily identified 

dec dec 
Resp 2 < ip* n+1 5L n A > = Res P2 < ip* +1 S(q n pl) A <fy>n,o > — 

= - < 0n+l,O<^n A { S Pl)^n,l > 
= I. 

At the point P 3 , y = 0(x~ N ~ 2 ), and 

= ^n,0X~ n + ^n,\X~ n ~ l H 

■071+1 = V'n+l.O 3 ' + + "0n+l,l X + ' ' ' 



5.49) 
5.50) 

5.51) 
5.52) 

5.53) 

5.54) 
5.55) 

5.56) 
5.57) 
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It follows that the residue is given by 

dec 

Resp 3 < if>* +1 5L n A5ip n > — = W n+lfi 5{q n pl) A <^n,l + V£+i,i%«Pn) A 5ip nfi ] (5.5 



We now make use of the equation (|5.51| ) to derive recursion relations between the coefficients 

C+1,0 = -Tpn,0QnPn, $1+1,1 = ^",0 ~ C,l9nPl- (5.59) 

They imply that 

V£+l,o?n = 0, V£+l,l9n = V»n,0?n (5-60) 

As a consequence, the first term on the right hand side of () simplifies to 

V&-l,0%nP») A ^n,l = <+l,Q<^n A p£<%»,l (5-61) 

Now recall that we introduced the coefficient j3 n by ^ n = f3 n q n . Comparing with the equation 
(), we obtain 

Pn = -p^V«,i (5-62) 

and the preceding term becomes 

C+i,o%nPD A # n ,i = -C+i,o<^n A % - VCn,o( 5 <?n A Spl)^ n>1 (5.63) 

On the other hand, p^ n ,o = 0, and the second term on the right hand side of () can be 
rewritten as 

C+1,1 S (lnPl) A 5ijj nfi = i)* n+l l q n 5p T n A 5^ nfi - ^* n+1>1 5q n A (^)^n,Q (5.64) 
Altogether, we obtain the following expression for the residue at P3 

Resp 3 < i/j* +l 5L n A 5^ > = II + III (5.65) 

or 

where the terms II and III are defined by 

II = -K + i,o(*?nA*pI)i,i + C + nfcA^)i i0 ] (5.66) 

HI = -{^* n+lfi 8q n A 5/3 n - ^* n+Xil qJPn A #n,o) (5.67) 

We claim that the term III can be simplified to 

HI = -5pl a 5q n (5.68) 

In fact, in view of the recursion relations ( |5.59| ) and the fact that ip n = (3 n q n , it can be 
rewritten as 

III = -ipn,oQn)(pl$qn A 5(3 n + 8pl A (5(3 n )q n + 5pl A (3 n 5q n ) (5.69) 

The first two terms on the right hand side cancel, since p^Qn = 0. As for the remaining term, 
we note that the normalization ip*ip n — 1 implies near P 3 

1 = ($V~ n + 0(x- n - l ))(f3 n q n x n + 0(x n - 1 )) = r nfl Pnq n + Oix" 1 ) (5.70) 
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from which it follows that ■?/>* Q f3 n q n = 1. The identity ( p.68| ) is established. 



Finally, it is readily seen that the remaining terms I and II combine into 

3 dx 
I + II = - £ Res Pa < r n+ Mn A 5p^ n > — (5.71) 

But the 1-form < ip^ +1 8q n A Sp^ip n > — is meromorphic on the space T, with poles only at 
the points P a above x = oo and Q a above x = 0. We can deform then contours and rewrite 
II+III as residues at Q a 

3 dx 
I + II = £ Res Qa < ij*n + Mn A 8p T n ^ n > — (5.72) 

a=l X 

At x = 0, we have Vn+i = V'n; anc ^ this expression is determined by the normalization 
condition (|4.31 ) on the matrix W. In terms of ip n , the normalization (|4.31|) can be restated 



as the normalization condition ip*{0)ipn = ^ as an identity between 3x3 matrices. Thus 
I + II = 3 < 5q n A 5p n >, and we obtain the final formula for the symplectic form u in terms 
of p n and q n 

N+l 

co = 2j2^A 5p n . (5.73) 

71=0 

6 Hamiltonian theory and Seiberg-Witten differential: 
The Odd Divisor Model 



The main difference between the even and the odd divisor spin models is in the parity of 
the divisor D of poles of the Bloch function ip n {Q). For the odd divisor spin model, D is 
essentially odd under the involution a : (x,y) —>■ (— x, y _1 ) in the following sense 

3 

[D] + [D a ] = K + 2j2Pa (6.1) 

a=l 

Here K is the canonical class, which is the divisor class of any meromorphic 1-forms on 
T. As in the case of the even divisor spin model, the relation ( |6.1|) is a consequence of the 
transformation of L(x) under a, which is in this case L(— x) = (L(x) _1 ) T . This implies that 
ipo{Q a ) an d ipo(Q)* are both dual Bloch functions for L(x), and thus 

r (Q) = MQ a )f(Q) (6.2) 

where f{Q) is a meromorphic function on T. But the zeroes of the dual Bloch function 
are exactly the poles of ipo{Q), while its poles are exactly the branch points of the surface 
T. Thus the preceding equation implies the following equation for divisor classes 

[branch points] - [D] = \D a ] (6.3) 

To determine the divisor of the branch points of T, we consider the differential dx, viewed 
as a meromorphic form on V. Since dx has a pole of order 2 at each P a , and a zero at 
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each branch point, we have [branch points] — 2 J2 a =i Pa = K, and the desired relation ( |6.1| ) 
follows. 

• We discuss briefly the direct and the inverse problems for the odd divisor spin system. 
Once the difference in parity of the divisor of poles of the Bloch functions is taken into 
account, the direct problem is treated in exactly the same way as before. As for the inverse 
problem, we need only a few minor modifications in expansions near the punctures Pi,P 3 , 
which we give (c.f. (|4.50| , |4.51|) now 



<j> n = x n e Bt (J2<f>n,k(Pi>~ k > Q -> Pi, (6-4) 
K = x~ n e xt ( £ ^ fc (P 3 )x-* ] , Q -> P 3 , (6.5) 

\k=0 ) 

They lead to minor modifications in the exact formulas for the Baker-Akhiezer function 



<f) n {t,Q) (cf. O): 



KAt, Q) = UQ) e{A{Q) - Za) m _ +nV + Zq) exp IJ^ndQo + m j (6.6) 

where dQ~ = dSli - dti^ and U~ = U - U a . 

We show that if the divisor D satisfies (|6.1| ), then the corresponding Baker-Akhiezer 



function satisfies the relation 

<j>* n (t,Q) = ( j ) T n (t,Qnf(Q), (6.7) 
where as before 0* are the rows of the matrix inverse to the matrix 

KM = <t>n,a(Pp) (6-8) 

Here the points P a (x) are the three preimages of x on T on different sheets. Of course, the 
matrix $„(x) does depend on the ordereing of sheets, but one can check that if for P a (x) 
we define 0* (P a ) as the corresponding row of the inverse matrix, then <p* n is well-defined. As 
before </>* has poles at all the branching points and zeroes at the points of the divisor D. 

To establish ( |6.7|) , we show that 

£ <P n , a (t, PjixWnA^ P l)f( P l) = <W (6-9) 

a 

Indeed, from (|6.4j ) and ( |6.5|) , it follows that the function 0„ )Q ,(t, Q)(pn,/3(Q cr ) f (Q) is holomor- 
phic everywhere except at the branching points (the poles and the essential singularities at 
the punctures P a over x = oo cancel each other; there are no poles at D and D° because 
f(Q) has zeros at these points). Therefore, the left hand side of the above equation is a 
holomorphic function of x (the poles at the branching points cancel upon the summation). 
Hence it is a constant, which can be found by taking x = 0. 

The uniqueness of <p n and the relation ( |6.7| ) implies as before that it satisfy the equation 



■>n+l 



L n (x)4> n , d t <p n = M n (x)<p n (6.10) 
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where L n and M n have the form ( |3.4| . |3-5| 



• We come now to the Hamiltonian structure of the odd divisor spin model. Recall that 
we had introduced the space Ai odd of spin chains. Solving the direct and inverse spectral 
problem as in the case of the even divisor spin model, we can identity Ai odd with the space 
of geometric data 

.M? dd ~{r,£>; [D} + [D^]=K + 2Y / P a } (6.11) 

a=l 

We can verify that the space on the right hand side is 2N + 1 dimensional, as it should be: 
there are N + 1 moduli parameters for the curve T, and N parameters for the antisymmetric 
divisor [D]. The same discussion as in §5.3 and §5.4 for the even divisor spin model shows 
that, in the present case, the only candidate for symplectic form is the form oom, restricted 
to the 2iV-dimensional phase space .M° dd defined by 

M otd = M odd n { Ujf = X | ( 02 ) 
The corresponding action and angle variables are now given by 

2N+1 

d\ (lh 0i= E / d < dd i 1<*<N (6.13) 
JA i j= i J 

where duf d and Af d are respectively a basis of odd holomorphic differentials and a basis 
of odd A-cycles. We have then as before 

N 

W(i) =J2 Sa 3 A( %" ( 6 - 14 ) 
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